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£h ! ABSTRACT 

£^ . It is well known that the density and anisotropy profile in the inner regions of a stellar 

system with positive phase-space distribution function are not fully independent. Here 
we study the interplay between density profile and orbital anisotropy at large radii in 
■ physically admissible (consistent) stellar systems. The analysis is carried out by using 

two-component (n-7,71) spherical self-consistent galaxy models, in which one density 
q , distribution follows a generalized 7 profile with external logarithmic slope n, and the 

5-H ■ other a standard 71 profile (with external slope 4). The two density components have 

c/5 \ different "core" radii, the orbital anisotropy is controlled with the Osipkov-Merritt 

d . recipe, and for simplicity we assume that the mass of the 71 component dominates 

the total potential everywhere. The necessary and sufficient conditions for phase-space 
consistency are determined analytically, also in presence of a dominant massive central 
black hole, and the analytical phase-space distribution function of (n-7,1) models, and 
' of rt-7 models with a central black hole, is derived for 7 = 0, f , 2. It is found that the 

, density slope in the external regions of a stellar system can play an important role in 

C*~) ■ determining the amount of admissible anisotropy: in particular, for fixed density slopes 

' in the central regions, systems with a steeper external density profile can support more 

radial anisotropy than externally flatter models. This is quantified by an inequality 
formally identical to the "cusp slope-central anisotropy" theorem (An & Evans 2006), 
OO ! relating at all radii (and not just at the center) the density logarithmic slope and the 

. anisotropy indicator in all Osipkov-Merritt systems. 
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1 INTRODUCTION Two general strategies can be used to construct a (self) 

consistent model, or check whether a proposed model is 



Observationally it is well established that elliptical galax- 
ies have dark matter halos, and also host central super- 



(self) consistent: they are commonly referred to as the "/- 

, , , to-p" and the "p-to-f" approaches, where / is the model 

massive black holes, these empirical tacts motivate the , .. n . „ m ■ r,™ n \ a 

. , , , TTT , , DF (e.g., see Bertin 2000, Bmney & Tremame 2008). An 

study ot multi-component dynamical models. When study- ir.if, , . ri<- ■ , , 

f example ol the first approach is the survey ot self-consistent 

mg dynamical models of stellar systems (single or multi- , ' , , , . , , , „ , . . 

v two-component galaxy models carried out by Bertin and 
component), the minimal requirement to be met by a phys- , , ,, , nl , , , 

' , co-workers, where the stellar and dark matter components 

ically acceptable model is the positivity of the phase-space . ., , , . . , , „ ,. „ 

. are described by two Dts ot the /oo family (e.g., Bertin & 
distribution function (DF of each distinct component. A ... , ^ n 1 

v ' . btiavefli 1984, Bertin et af. 1992); other welf known exam- 

model satisfying this essential requirement (which is much , ,, T ,. / lr> ^^\ 11 1 ^ ? 1 i /t-. j.- 

. .. . . pies are the King (1966) models and the j v models (Bertin 

weaker than stability, but stronger than the fact that the „ m .. „..„, TT - , , , ,, .. , , 

„ & Irenti 2003). Unfortunately, the resulting spatial densi- 

Jeans equations have a physically acceptable solution) is .. , , . , , , . . . , „ . 

' ties obtained by solving the associated Poisson equation are 

called consistent moreover, when the total gravitational po- . , . „ . , , 

j. t . , . j + J,, j u ii. * j. 1 j •* ci-tu u m S enera i n °t expressible m terms of simpfe or even known 

functions, and so only numerical investigations are usually 

feasible. In the "p-to-f" approach the density distribution 

is given, and specific assumptions about the model internal 
gravitating system. , . , . . , . . „ , r 

dynamics are made; m special cases inversion formulae from 

the density to the DF can be obtained, usually in integral 

* „ „, , y j. I" - t~\ t — 1 1 "i /— t ■ form or series expansion (see, e.g., Fricke 1952, Lynden-Bell 

Current address: Max-Planck-lnstitut lur bx. Pnysik, Giessen- v 

bachstrafie, D-85741 Garching, Germany 1962 > 0si P kov 1979 > Merritt 1985 > hereafter OM; Dejonghe 



tential is determined by the total density profile through 
the Poisson equation, the model is called self-consistent. 
In other words, we call self-consistent a consistent self- 
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1986, 1987; Cuddeford 1991; Hunter k Qian 1993, Ciotti 
& Bertin 2005). In particular, in order to recover the DF 
of spherical models with orbital anisotropy, the OM tech- 
nique has been developed from the Eddington (1916) inver- 
sion for isotropic systems, and widely used to study one and 
two-component models (see, e.g., Ciotti & Pellegrini 1992, 
hereafter CP92; Hiotelis 1994; Carollo et al. 1995; Ciotti & 
Lanzoni 1997; Ciotti 1996, 1999, hereafter C96, C99; Baes 
& Dejonghe 2004; Buyle et al. 2007). We remark that the 
OM parameterization is not necessarily the best description 
of real systems, however its simplicity and the fact that it 
captures the main features of models of galaxy formation, 
that are generally found nearly isotropic at the center and 
increasingly radially anisotropic in the outer envelope (e.g., 
van Albada 1982; Trenti, Bertin & van Albada 2005; Nipoti, 
Londrillo & Ciotti 2006; Binney & Tremaine 2008. But see 
Cuesta et al. 2008 and references therein) , make it the natu- 
ral choice for investigations as that presented in this paper. 

In many cases, the difficulties inherent in the operation 
of recovering analytically the DF prevent a simple consis- 
tency analysis, and phase-space positivity must be inves- 
tigated by numerical inspection of the inversion integral. 
In these cases the reasons underlying consistency or incon- 
sistency of a proposed model tend to be obscured by the 
numerical nature of the solution. Fortunately, informations 
about consistency of multi-component OM systems can be 
obtained without recovering their DF, following the proce- 
dure described in CP92. This method uses the radial density 
profile of each component and the total potential of the sys- 
tem, and gives necessary and sufficient conditions for (self) 
consistency. Moreover, since only spatial differentiation and 
inequality checks are required, this method is best suited for 
analytical investigations. For example, C96 and C99 applied 
the CP92 technique to the general family of two-component, 
spherically symmetric and radially anisotropic (71,72) mod- 
els. This family is made of the superposition of two 7 mod- 
els (Dehnen 1993, Tremaine et al. 1994) with different total 
masses, scale-lengths, inner density slopes, and OM radi- 
ally anisotropic velocity dispersions. The possibility to in- 
vestigate the combined effects of radial anisotropy and inner 
density slope on multi-component systems made the study 
of (71, 72) models interesting, as it is well known that the in- 
ner density profile sets important constraints on the amount 
of admissible radial anisotropy (e.g., Richstone & Tremaine 
1984), and indeed in C96 and C99 analytical limitations on 
anisotropy as a function of the density slopes 71 and 72 were 
obtained. These models clarified the reasons behind the nu- 
merical findings of CP92, i.e. the difficulty of consistently 
superimposing a centrally peaked distribution such as the 
de Vaucouleurs (1948) profile to a centrally flat one, such as 
the King (1972) or quasi-isothermal density profile (even in 
the isotropic case). In fact, it was shown that the DF of the 
71 component in isotropic (71,72) models is nowhere nega- 
tive, independently of the mass and concentration of the 72 
component, whenever < 72 < 71 and 1 < 71 < 3. On the 
contrary, a 71 = component becomes inconsistent when 
adding a 72 = 1 halo with a small core radius. Thus, in 
two-component isotropic models, the component with the 
steeper central density distribution is usually the most ro- 
bust against inconsistency. More recently, the importance of 
the central density slope in limiting the amount of possible 
radial anisotropy has also been quantified with the so called 



"cusp slope-central anisotropy" theorem (An & Evans 2006, 
hereafter AE06; see also eq. [28] in de Bruijne et al. 1996). 

The previous investigations left however unexplored the 
importance of the external density slope in determining the 
model consistency. In fact, the phase-space density cannot 
be identified, in general, with any specific spatial position in 
the system, as (for example) stars of a given energy can span 
a large range of radial position^] (systems made of circular 
orbits are an obvious exception) . Therefore also the external 
regions of a density distribution can be important in limiting 
the maximum allowable anisotropy, but the (71,72) models 
are of no help in the study of this issue, because the exter- 
nal density profiles of both components all decrease as r -4 . 
For these reasons here we focus on the phase-space prop- 
erties of n-7 models, i.e. models similar to the standard 7 
models in the inner regions, but with a density profile pro- 
portional to r~" (instead of r~ 4 ) in the external regions; re- 
markably, several properties of n-7 models can be obtained 
from those of 7 models by differentiation with respect to 
their scale-length. In this notation, the 4-7 models coincide 
with the standard 7 models. We also study the larger class 
of (n-7,71) models, i.e. two-component systems in which a 
71 halo is added to a n-7 component. Thus, here we further 
explore the trends emerged in CP92, C96, and C99, deter- 
mining the limits imposed by phase-space consistency on 
the parameters describing (n-7,71) models, and n-7 models 
with a central BH [hereafter (n-7,BH) models], with partic- 
ular focus on the effects of the external slope parameter n. In 
specific cases (that we call halo-dominated models), the cal- 
culations are performed under the assumption that the mass 
of the halo component (or of the central BH) is dominant 
over the mass of the visible one. This assumption is mainly 
motivated by mathematical simplicity (see also Sect. 4), al- 
though this is not the only reason. In fact, for any given 
two-component model, it is expected that the DF proper- 
ties are bracketed by those of the one component model and 
by those of the halo-dominated model (corresponding to the 
formal case of an infinite halo mass). Of course, while the 
case of dark matter dominated systems can be considered 
a viable representation of some real astrophysical systems, 
the case of a dominant BH is less natural, and it just gives 
the strongest possible limitations for consistency of systems 
with a central BH. 

The paper is organized as follows. In Section 2 we briefly 
review the technique developed in CP92, and we prove that 
the necessary condition for consistency derived in CP92 for 
OM systems can be rewritten exactly as the AE06 "cusp 
slope-central anisotropy" theorem, holding however at all 
radii and not just at the center. In Section 3 the one and two- 



1 Incidentally, this implies that the use of p/cr 3 as a proxy for 
phase— space density, where a is the local value of the velocity dis- 
persion, has no assignable meaning without an appropriate discus- 
sion. For example, in power— law isotropic systems with p ex r~ 7 , 

■y-6 

because the functions pi a cc r 2 (for 7 > 1) and phase-space 

6-7 

density f(£) oc £ 2 (~<- 2 ) (for 7 > 2, where £ = —E is the so— called 
relative energy, see Sect. 4 and Baes et al. 2005) are both power 
laws with respect to their arguments, the exponents are related 
in a simple way. However, the converse statement is not true: for 
example, in the Plummer (1911) sphere /(£) oc £ 7 / 2 is a power 
law at all energies, but p/cr 3 is not a power law of radius. 
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component n-7 models are introduced, and the necessary 
and sufficient conditions imposed on the model parameters 
by phase-space consistency are derived for different values 
of the logarithmic density slope n. In Section 4 the DF of 
the n-7 component of halo-dominated (n-7,1) models with 
7 = 0, 1,2, and of similar models with a dominant central 
BH, are derived explicitly for arbitrary (but integer) values 
of n, and the true boundaries of the consistency region in 
the parameter space are obtained. The main results of the 
investigation are summarized in Section 5. Finally, in the 
Appendix an easy method to solve analytically the Jeans 
equations in the general OM case for the wider class of (ni- 
71,^-2-72) models is presented. 



2 THE CONSISTENCY OF 

MULTI-COMPONENT OM SYSTEMS 

A stellar system made of the sum of the density compo- 
nents pit is called consistent if each DF (f k ) is non-negative 
over the whole accessible phase-space; a consistent self- 
gravitating system is called self-consistent. The technique 
developed in CP92 permits us to check whether the DF 
of a multi-component spherical system, where the orbital 
anisotropy of each component is modeled according to the 
OM parameterization, is positive; in practice, only infor- 
mation about the radial trend of each density component 
and of the total integrated mass are used. In the OM for- 
mulation the DF of each component is obtained assuming 
/k = /k(Qk), with 



= £ 



J 2 



(1) 



where £ and J are respectively the relative energy and the 
angular momentum modulus per unit mass, r a k is the so- 
called anisotropy radius, and /k = for Q k < (e.g. see 
Binney & Tremaine 2008). The velocity dispersion tensor 
associated with eq. ([T} is characterized by radial anisotropy 
increasing with the radius r, while in the limit r a k — > 00 the 
system becomes globally isotropic. As well known, the DF 
of the density component pk is given by 



/k(Qk) = 



_L d_ 



1 



V8tv 2 



VOk - *T 

d 2 g k d^T 



where the augmented density is 



£>k(r) 



1 + 



Pk(r), 



(2) 



(3) 



^ T ( r ) ~ X^fe ^k{r) is the total relative potential, < Qk < 
^I/t(0), and in the integrals above it is assumed that the 
radius is eliminated from g k in favour of $t- It can be proved 
that the second equivalence in eq. ([5} holds for untruncated 
systems with finite total mass, as the models discussed here. 



2.1 Necessary and sufficient conditions for 
consistency 

Quite obviously, only a very limited number of density pro- 
files among those expressible in analytic form admit an 



explicit DF, so that the study of phase-space consistency 
would appear restricted to such rare cases when conducted 
analytically, while all the remaining cases should be investi- 
gated numerically. Fortunately this is not true, as the CP92 
technique is based on the verification (numerical or analyt- 
ical) of the following 

Theorem: A necessary condition (NC) for the non- 
negativity of /k is 



dg k {r) 
dr 



< 0, < r < 00. 



(4) 



If the NC is satisfied, a strong sufficient condition (SSC) for 
the non-negativity of /k is 



d_ 

dr 



dg k (r) rV*T(r) 



dr 



Mt(t-) 



> 0, < r < 00. 



(•>) 



Finally, a weak sufficient condition (WSC) for the non neg- 
ativity of /k is 



d_ 

dr 



dg k (r) 



dr M T (r) 



> 0, < r < 00. 



(6) 



Proof: See CP92, C96, C99. 

The first consideration that follows from the conditions 
above is that the violation of the NC is connected only to 
the radial behavior of the augmented density g k , and so this 
condition applies independently of any other density com- 
ponent of the model. Obviously, a model failing the NC is 
certainly inconsistent, while a model satisfying the NC may 
be inconsistent, i.e. the f k may be negative even for values 
of model parameters allowed by the NC. The second con- 
sideration is that a model satisfying the WSC (or the more 
restrictive SSC) is certainly consistent, while a model failing 
the WSC (SSC) may be consistent. Therefore the consistency 
of an OM model satisfying the NC and failing the WSC (or 
the SSC) can be proved only by direct inspection of its DF. 



2.1.1 A density slope-OM anisotropy inequality 

The NC can be recast into a simple inequality between the 
value of the density slope 



7(r) 



dlnp 



din r 

and the value of the orbital anisotropy indicator 



2a? 



r + r 2 



(7) 



(8) 



that must hold at each radius for consistent OM systems. In 
the expression above a 2 and a 2 are the radial and tangential 
components of the velocity dispersion tensor of the system 
(e.g. Binney & Tremaine 2008), and in the following proof we 
restrict to natural systems, i.e. systems with monotonically 
decreasing density profile, so that 7(7*) > 0. The proof is 
trivial: in fact, it suffices to express the NC in terms of the 
logarithmic density slope as 



2 
7(r) 



<0, 



r 2 + r 2 

and from eq. (J8]l one obtains the necessary condition 
7 (r) > 2/3(r), 



(9) 



(10) 
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which must be verified by each OM component of a con- 
sistent multi-component system. Curiously, The condition 
above is formally identical to the inequality appearing in the 
"cusp slope-central anisotropy" theorem derived in AE06. 
This latter theorem was proved by using constant anisotropy 
systems (i.e., /3(r) = /3, and in this case as well the inequal- 
ity above holds at all radii in consistent models. See eq. [10] 
in AE06), and then it was convincingly argued that the in- 
equality it holds asymptotically for the central regions of 
spherical systems with generic anisotropy distribution. 

In the specific case of n-7 density distributions (see 
eq. [IS]), the logarithmic density slope 

and therefore eq. (|10|l shows that in general the NC is satis- 
fied near the center of OM systems (where /3 = and j(r) = 
7), and at large radii (where (3 — 1 and j(r) — n > 3), while 
critical behaviors may be expected at intermediate radii (see 
Sect. 4). 

2.2 Classification of phase-space inconsistency as 
a function of the anisotropy radius 

Interestingly (but not unexpectedly) the particular func- 
tional form of the augmented density characteristic of OM 
parametrization, limits the possible manifestations of phase- 
space inconsistency to few general cases, that can be illus- 
trated as follows. From eqs. ©-© it is apparent that the 
DF of each density component can be written as 

U(Q) 



HQ) = fi(Q) 



(12) 



where f\ is the DF of the considered density component in 
the isotropic case (for simplicity, from now on we avoid the 
use of the subscript k). Let A+ be the subset of phase-space 
defined by the property that fx is positive VQ € A+. Then, 
from eq. (|12p a first condition to be satisfied for consistency 



> rl 



\ 



0, sup 



HQ) 



HQ) 



(13) 



Obviously, when /j > over all the accessible phase-space 
(the common situation for realistic density distributions), 
inequality (113[) is also the only condition to be satisfied for 
the model consistency, and there is at most a lower bound 
for the anisotropy radius. For example this is the case for 
the families of one-component anisotropic 7 models (C99) 
and Sersic (1968) models (Ciotti 1991). 

When the complement of A+ is not empty, i.e. /; < 
over some region A- of the accessible phase-space, a second 
inequality again derived from eq. (|12p must be verified for 
consistency: 



<r2 



inf 



.HQ) 



\MQ)\ 



(14) 



Q6A_ 



Therefore, if there exists some Q £ A- for which / a < 0, 
then the proposed model is inconsistent^. If / a is positive 



2 In C99 and Ciotti (2000) it is erroneously stated that the model 
is inconsistent if / a < VQ € A—. All the results presented 
therein are however correct. 



over A-, from conditions (|13[) - (|14| ) it follows that r~ < r a < 
r£ for OM consistency, so that if < r& the proposed 
model is inconsistent. Note that formally identical consid- 
erations hold when discussing the inequalities (|4]) and ©, 
because from eq. (|3]) it follows that the NC and WSC can be 
written in the same functional form of eq. (|12jl : of course, 
in these cases the sets A+ and A- are to be intended as 
subsets of the radial range < r < 00. 

In the following Section, after presenting the one and 
two-component n-7 models, the necessary and sufficient 
conditions for model consistency will be derived, also for 
the case of n-7 models with a massive central black hole. 



3 THE ONE AND TWO-COMPONENT 
n-7 MODELS 

The one-component n-7 models are a natural generaliza- 
tion of the widely explored family of the so-called 7 models 
(Dehnen 1993, Tremaine et al. 1994), whose density, cumu- 
lative mass within r, and relative potential are given by 

3-7 Mr c 



p(r) = 



47T r~'(r c +r) 4 "T' 
3-7 

MM = M x ( — — 

w \r c + r 



*(r) = 



GM 



(2-7) 
r + r c 



2- 7 



(15) 



(16) 



(17) 



hi 



, (7 = 2), 



respectively. 

In the n-7 models the logarithmic density slope for 
r > r c is not 4, but it is a free parameter n > 3 (to en- 
sure that their total mass is finite), so that these density 
profiles belong to the family considered by Zhao (1996). In 
the following, in order to exploit a useful analytical property 
of the n-7 density profiles, we will often assume n restricted 
to integer values > 4. In fact, the generic density profile of 
a n-7 model is given by 



Pn(r) : 



Mr" 



47rB(3 -7,n - 3)rT(r c + r) n ~~i 



(18) 



r?~ 3 (-l)" rf"~ 4 p(r) 
T(n-3) dr?- 4 r c : 



where M is the total mass of the density distribution, 
B(s,y) = T(x)r(y)/T(x + y) and T(x) are the complete 
Euler beta and gamma functions, respectively, and the first 
expression holds for any real number n > 3. The second ex- 
pression, based on repeated differentiation of eq. (|15|) with 
respect to r c , holds instead for integer n > 4. Of course, for 
n = 4 the standard 7 model density profile is recovered. The 
radial behaviour of p n (r) is shown in Fig. [1] (top panel) for 
n-1 models with increasing n. In the bottom panel we show 
the corresponding logarithmic density slopes, calculated ac- 
cordingly to eq. ([TT}. In particular, it is apparent how the 
inner slope is 7 (for r<"/r c /n), and the external is n (for 
r>r c ). As a consequence, while the density profile outside 
r c becomes more and more steep at increasing n, the inner- 
most region where the density slope is 7 contracts near the 
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L °glo( r / r c) 



Figure 1. The normalized density profile p n of the n-1 model 
(top panel), and its logarithmic slope (eq. 1111 , bottom panel), 
for n = 4,5, and 6. 



center. As anticipated, the possibility to write the density 
p n as a derivative with respect to the scale-length r c pro- 
vides an easy way to determine analytical properties of the 
n-7 models for n integer. Indeed, since the density enters 
linearly in the integrals of the total mass inside the radius r 
and of the potential, from eq. {TSJ it follows that 



M„(r) = 



*n(r) = 



r c ra ~ 3 (-l)" rf"~ 4 M(r) 
r(n-3) dr?" 4 r c ' 



(19) 



(20) 



r(n-3) dr?" 4 r c 

Instead, when n is not at integer, M n and ty n are in general 
given by hypergeometric 2-P1 functions. Expressions similar 
to eqs. (|19|l - (|20p hold for any quantity that can be written as 
a linear functional of p n , so that the surface density profile of 
n-7 models with n integer can be obtained by repeated dif- 
ferentiation of the surface density profile of 7 models (when 
analytically available, e.g. see Binney & Merrifield 1998). 
This property will be exploited in Section 4 to obtain the 
explicit DF of halo dominated (n-7, 71) models; moreover, in 
Appendix B we show how quadratic functionals of p n (such 
as the gravitational energy and the velocity dispersions) can 
be also evaluated by using repeated differentiation with re- 
spect to r c . 

We are now in the position to introduce the two- 
component models used in this work. The most general 
family of two-component n-7 models, i.e. the (711-71,712-72) 
models, is made of the superposition of two n-7 models with 



different total masses, scale-lengths, internal and external 
density slopes, and finally two different anisotropy radii. For 
simplicity here we restrict to the case of (71-7,71) models, 
where the "halo" density distribution is a standard 71 model: 
we note however that some of the presented results can be 
generalized without much effort to the family of (711-71,712- 
72) models (see Appendix B). In the following, the total 
mass Mi and the characteristic scale-length r\ of the 71 halo 
are adopted as normalization constants, so that the physical 
scales for density and potential are given by pN = Mi/rf 
and = GMi/ri, while we define s = r/ri, £ = r c /r\, 
p = M/Mi. With this choice, the (71-7,71) models are struc- 
turally determined by fixing the four independent parame- 
ters (Mi, n, fi, £), with the obvious conditions p > and 
£ > 0. We conclude this introductory discussion by noticing 
that, for reasons that will become apparent in Section 3.2 
and 4, the present normalization differs from that adopted 
in C99, where the normalization mass and scale-length were 
those of the first component. 



3.1 The necessary and sufficient conditions for 
one— component n-7 models 

Before discussing the case of the two-component models, 
we consider the NC for anisotropic n-7 models, in order to 
determine analytically a critical value for the anisotropy ra- 
dius such that a higher degree of radial OM anisotropy (i.e., 
a smaller r a ) would produce a negative DF for some permit- 
ted value of Q. We recall that the obtained anisotropy limit 
holds also when a second component is added (see Section 
2), so that the present discussion is fully general. Moreover, 
we note that as the NC involves only the density distribution 
under scrutiny, we can use the first expression in eq. (|18p . 
and the following results hold for any n > 3, not necessarily 
limited to integers. In the following the unit mass and unit 
length are the total mass M and the scale-length r c of the 
n-"f model, with s a = r a /r c . 

In C99 the analytical expression for the critical s a (7) 
was obtained for the whole family of 7-models, and here we 
derive its generalization s a (7, n). In fact, eq. (Al) shows that 
for 2 < 7 < 3 and n > 3 the NC is satisfied for s a > 0, and 
the result of C99 is now obtained as the particular case for 
n = 4. In other words, the NC leaves open the possibility of 
making n-7 models with 7 > 2 using radial orbits only. In 
the range < 7 < 2 the NC requires instead that 



> sm 



7 - (n - 2)s M 



7 + nsM 



(21) 



where the explicit expression of 7) is given by eq. (A2), 

and the inequality above reduces to eq. (13) of C99 when 
n = 4. The NC then proves that n-7 models with < 7 < 
2 cannot be made of radial orbits only, independently of 
the value of the external slope n, and of the presence of 
any possible second component. This result extends the list 
of cases for which it has been proved that a density cusp 
shallower than r~ 2 cannot be supported by radial orbits 
only (Richstone & Tremaine 1984, CP92, C99, AE06). In 
Fig. [2] the lower bound for the anisotropy radius given in 
eq. (|21|l is shown with the solid lines as a function of 7 for 
n = 4, 5, 6: n-7 models (one or multi-component) with the 
pair (s a ,7) in the nearly triangular region under the solid 
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Figure 2. Consistency limits on s a = r a /r c as a function of 7 
in one— component n-7 models, for n = 4,5, 6. The NC limit is 
represented by the solid lines: all models with the pair (s a ,7) in 
the triangular region below them are inconsistent. The dashed 
lines mark the WSC limit: all points above these lines correspond 
to consistent one— component n-7 models. For 7 = 0, 1, 2, the filled 
circles joined by the dotted line represent the more accurate limits 
obtained from the SSC. For the black-hole dominated (n-7,BH) 
models the long-dashed lines (interrupted for 7 < 1, see Section 
3.2), are the lower limit on s a given by the WSC, while the solid 
squares are the true limits derived from the DF. 

curves are inconsistent. In particular, for fixed 7 an increase 
of n produces a decrease of the minimum s a , i.e. we have 
here a first indication that a steepening of the density profile 
in the external regions of a system with fixed inner density 
slope can be effective in increasing the maximum amount of 
sustainable radial anisotropy. This behaviour is quantified 
by substitution in eq. (|21|l of the asymptotic expansion for 
n — > 00 of sm(Vi,7) in eq. (A2): 

/ x l-2 7 + VI +47 _ 2 

SM(n,7) = — hO(n ). (22) 

As often happens in asymptotic analysis, even if the expan- 
sion above holds in principle only for very large values of n, 
the substitution of (|22p in (|21|) leads to percentual errors on 
s a less than 22%, 14%, and 10% for n — 4,5,6 respectively 
(for the inner density slope 7=1). 

We now move to discuss the WSC for one-component 
n-7 models: the obtained s a (fi, 7) will mark a lower limit 
above which consistency (for the considered n and 7) is 
guaranteed. Unfortunately, the WSC cannot be explored 
algebraically in the general case, because the resulting in- 
equality (that for simplicity we do not report here) involves 
the solution of an equation of fifth degree for n = 5, and 
the degree increases for increasing n. For n — 4, instead, it 
is possible to treat the WSC analytically, since it reduces 



to the discussion of a cubic equation (see C99). Of course, 
the critical values for the anisotropy radius can be easily ob- 
tained solving numerically inequality © in any specific case 
of interest, and the results are shown in Fig. [2]for n — 4,5,6 
with dashed lines: all one-component n-7 models with the 
pair (s a ,7) in the region above the dashed lines are consis- 
tent. Note how for increasing values of the external density 
slope n more and more radially anisotropic orbital distri- 
butions can be supported, thus confirming the indications 
provided by the NC. 

Values of s a nearer to the limits set by the DF are 
obtained by using the SSC. Inequality (JS| evaluated for a 
generic pair (7,71) results in a transcendental equation that 
must be studied numerically (as already done for n = 4 in 
C99), and the black dots joined by the dotted lines in Fig. [5] 
represent the critical lower values of s a for one-component 
7 models with 7 = 0,1,2, and n = 4,5,6. As expected, 
the dotted lines are contained between the solid (NC) and 
the dashed (WSC) lines in each panel, and again they shift 
downward for increasing n. 

Thus, from this preliminary investigation of one- 
component n-7 models we conclude that for fixed 7 an in- 
crease of n corresponds to a decrease of the minimum ad- 
missible value of s a , i.e. steeper density distributions in the 
external regions can support more radial anisotropy. 

3.2 Sufficient conditions for halo— dominated 
(n-7,71) models 

In order to proceed further with the preliminary discussion, 
we now apply the WSC to the n-7 component of a (n-7,71) 
model, extending to this class of systems some of the re- 
sults obtained in C96 and C99 for two-component (71,72) 
models. In particular, the analytical study in C99 was re- 
stricted to some representative (71,72) models, namely a) 
isotropic two-component systems with inner slopes in the 
ranges 1 < 71 < 3, and < 72 < 71 (i.e. the 72 compo- 
nent is shallower in the central regions); 6) isotropic two- 
component (0,1) systems (i.e. the 72 component is steeper 
in the central regions); and finally c) anisotropic 71 profiles 
with 1 < 71 < 3 in the gravitational field of a dominant 
central black hole. Here, in order to obtain analytical re- 
sults for the more general (n-7,71) models, we assume that 
the 71 component (the "halo") dominates everywhere the 
gravitational field. Under this simplifying assumption, the 
following three results, corresponding to the points a), b), 
and c) above, will be proved analytically: 

(i) In the case of halo-dominated isotropic (n-7,71) mod- 
els, with 1 < 7 < 3, < 71 < 7, and n > 3, the centrally 
more peaked n-7 component is consistent, independently of 
the value £ = r c /?T of its concentration relative to the cen- 
trally flatter 71 halo. In the case of anisotropic (n-7,71) mod- 
els, the determination of a minimum anisotropy radius for 
consistency as a function of n,7,7i,£ reduces to the solu- 
tion of an algebraic equation of sixth degree (for generic 
n), which is solved numerically. In the particular case of 
halo-dominated (n-2,1) models, the application of the WSC 
shows that for £ < (n— 1)/2 these models can be consistently 
assembled using radial orbits only. 

(ii) In the case of halo-dominated isotropic (n-0,1) mod- 
els, the WSC shows that for £ < (n + l)/2 the n-0 density 
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Figure 3. Minimum value of the anisotropy radius (normalized 
to the scale-length of the n-7 component, i.e. s a = r a /r c ) as a 
function of the relative concentration £, for consistency of halo 
dominated (n-7,1) models. The external density slope n increases 
as n = 4, 5, and 6 from left to right panels, while the internal 
density slope as 7 = 0, 1, and 2 from top to bottom panels. Solid 
lines represent the limits derived from the WSC, while dashed 
lines are the Sa derived from the DF. Note that for (n-0,1) models 
also the Sj limits appear for £ > £c (see Sections 3.2 and 4.3 for 
details). 



distribution is consistent. For broader n-0 density distribu- 
tions, the models can be consistent only in presence of some 
amount of radial anisotropy. 

(iii) In the case of anisotropic n-7 models with a domi- 
nant black hole at their center it is possible to determine 
analytically a lower limit r a (n,7) for consistency, and this 
limit decreases for increasing n. 

The proof of the first result, which is an extension of the 
study mentioned in point a), is conceptually straightforward 
but algebraically cumbersome, and only the main steps are 
reported in Appendix A (eq. [A3]). In particular, as this 
result holds also for n > 3, while the external logarithmic 
slope of the 71 halo is 4, it means that a centrally steep 
density profile can be consistent in the gravitational field of 
a more massive (but centrally shallower) component, even 
when its external regions are less peaked. The situation is 
illustrated for some representative cases in Fig. [3] (bottom 
and middle panels), where all points above the solid lines 
correspond to consistent halo-dominated (71-1,1) and (n-2,1) 
models, and apparently the isotropic limit (s a = 00) is al- 
lowed for any choice of 5; the solid lines are obtained by solv- 
ing numerically the corresponding WSC. A few additional 
trends are apparent. First, for a given central logarithmic 
density slope 7 and for a given relative concentration 5, an 
increase of n corresponds to a better and better ability to 



sustain radial anisotropy, and thus the trend found in one- 
component models is confirmed also in the two-component 
case. Second, at fixed n and 7, the minimum anisotropy ra- 
dius for consistency increases for increasing 5- This trend 
was already found in C96 and C99: in practice, in a fixed 
potential broader density distributions are less and less able 
to sustain radial anisotropy. Third, at fixed external slope 
n and relative concentration 5, more centrally peaked sys- 
tems are better able to support radial anisotropy. An addi- 
tional comment concerns the specific case of (n-2,1) mod- 
els (Fig. [3] bottom panels). In fact, it is apparent how in 
presence of a dominant 71 = 1 halo the WSC limits flat- 
ten to zero for relative concentration less than some critical 
value 5c, and accordingly the model may be purely radi- 
ally anisotropic. Remarkably, it is easy to show that the n-2 
component can support purely radial orbits for relative con- 
centrations 5 < 5c = (n — l)/2 (see eq. [A4]), which are 
exactly the critical points in Fig. [3] Of course, this is just 
a sufficient condition, so we expect that the existence of a 
larger critical concentration for consistency in case of purely 
radial orbits will be revealed by direct inspection of the DF. 
A final comment is due. In fact, the adopted range of values 
for the internal density slopes 7 > 71 imply that the halo 
density distribution is centrally less peaked than the n-7 
component. Therefore, for large but finite total halo mass, 
the integrated mass of the halo is subdominant with respect 
to the stellar one for vanishing r, and the assumption of a 
dominant halo breaks down at the very center. This means 
that in an asymptotic sense, the analysis presented in the 
previous Section should be applied at the center of these 
finite-halo mass models. 

Result (ii) extends the study mentioned in point 6) 
above, and it is proved in eq. (A5). The obtained limita- 
tion 5 < 5c = (n + l)/2, shown by the solid line in Fig. [4] 
suggests that the concentration of the 7 = component 
must "adapt" to the dominant 71 = 1 halo for phase-space 
consistency. Note that 5c increases for increasing n, i.e. a 
steeper external density profile can (partially) compensate 
for the effect of a central shallower density distribution. The 
appearence of the concentration limit 5c in the isotropic case 
manifests in the top panels of Fig. El In fact, as discussed in 
Section 2, in systems with positive isotropic functions (WSC 
or DF) only a lower limit for s a exists, as in the case of halo 
dominated (n-1,1) and (n-2,1) models. In the present case, 
instead, when for increasing 5 the critical value 5c = 2.5 
is crossed (in the n = 4 panel), the additional condition 
s a < corresponding to the radial domain A- appears as 
the vertical asymptote of the solid line, and the parameter 
space for consistency progressively reduces, shrinking to zero 
for the value of 5 where the two solid lines (obtained numer- 
ically) cross. A similar configuration repeats in the panels 
representing the n = 5 and n = 6 cases, where the consis- 
tency limits move, as expected, downward and rightward. Of 
course, these limits represent only a sufficient condition for 
the consistency, and when considering the DF in the limit 
fi — > (see Section 4) we expect to determine larger criti- 
cal values of 5c, and larger consistency regions in the (5,Sa) 
plane. 

The result concerning the models with central black 
hole, related to point c) above, is proved as follows. As well 
known, for n — > the potential of the 71 = 1 model becomes 
that of a point mass, and from the first result proven in this 
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Section it follows that a black hole of any mass Mbh can 
be consistently added at the center of a globally isotropic 
n-7 model, when 1 < 7 < 3 and n > 3. For this reason 
in the following we restrict to this range of slopes 7, and 
accordingly the long-dashed curves in Fig. [2] interrupt for 
7 < 1. Having reduced the study to models with a positive 
isotropic DF, we only have to estimate a lower limit of r a . In 
Appendix A the WSC is applied to the anisotropic (n-7,BH) 
models with 1 < 7 < 3; it can be discussed analytically in 
the special case of a dominant black hole, i.e. assuming in 
eq. (O Mr = Mbh (and so * T = GM BH /r). As shown in 
eq. (A6), under these assumptions 



Sa > Sm 



(3 - 7) (7 - 2) + s M (n - 2)[(6 - 2 7 ) - (ra - 3)s M ] 



n(n - l)sjj + 2n(7 - l)s M + 7(7 - 1) 



(23) 



where sm = sm (n, 7) is obtained by solving an algebraic 
equation of fourth degree: the obtained limits are shown in 
Fig. [2] with the long-dashed lines, for n = 4, 5, 6. For n — 4 
equation above coincide with that given in eq. (18) of C99, 
while asymptotic analysis proves that for n — > 00 



, , S M o(7) , r., -2n 

s M (n,j) = h 0(n ) 



(24) 



where smo(7) is the larger real root of eq. (A9). We found 
that the substitution of eq. (|24[) in eq. (|23p leads to estimates 
of the minimum s a discrepant from the values obtained by 
the full formula by less than 50%, 33%, and 23% for n — 
4, 5, 6 respectively (for the inner density slope 7 = 1). As in 
the other cases, an increase of the external density slope n 
makes it possible to sustain more radial orbits. 



4 THE DF OF HALO-DOMINATED (n-7,71) 
MODELS 

It should be clear that the DF of two-component (n-7,71) 
models cannot be obtained analytically, except for very spe- 
cial combinations of the values of 71,7, and 71, as for ex- 
ample the (1,1) models discussed in C96 (that would be re- 
ferred as (4-1,1) models in the present notation). However, 
the search is not hopeless. In fact, from C99 it is known 
that the halo dominated DF of two-component (1,0) and 
(0,1) galaxy models can be obtained in terms of elemen- 
tary functions, and a simple argument shows that this is 
also the case of (2,1) and (2,0) models. This fact, the lin- 
earity of the OM inversion formula with respect to p in the 
non self-gravitating case (see eq. [2]), and finally the for- 
mula {18} for n integer, prove that also the DF of the halo 
dominated (n-7,1) and (ra-7,0) models can be expressed in 
terms of elementary functions, for any integer n > 4 and 
7 = 0, 1, 2. Note that, albeit their special nature, the study 
of halo-dominated models is useful because the formulae 
- expressible using elementary functions - can be studied 
very easily, making apparent the effects of the relative dis- 
tribution of the investigated component and of the halo. In 
particular, in this Section we will determine the DF of the 
n-7 component of halo dominated (n-7,1) models, and the 
exact phase-space constraints will be derived and compared 
with those obtained using the NC, WSC, and SSC in Sec- 
tion 3. We restrict to the 71 = 1 case because the Hernquist 



(1990) potential is the simplest in the class of 7 models; fur- 
thermore, this choice allows us to investigate the consistency 
of models with central density slopes flatter (7 = 0), equal 
(7 = 1), and steeper (7 = 2) than that of the halo. As a 
consequence, the results obtained for (n-7,1) models should 
be representative of the whole situation. 

With the normalization scales introduced at the begin- 
ning of Section 3, the relative potential of the Hernquist halo 
needed for the recovery of the DF is 

*(r) = = *n*(s), (25) 

1 + s 

while from eq. (I18|l 



Pn{r) = ON 



Mr- 3 (3-7)(-l) n d"- 4 ($ + s) 



7 — 4 



(26) 



47rr(n - 3)st d£™" 4 

As for the density and the potential, also for the DF it is 
useful to work with dimensionless functions, and we define 
/ = /n/(<3) with / N = Pn*n 3/2 > and < = Q/*n < 1: 
note that, at variance with C99, the normalization quantities 
used here are those of the halo, consistently with the halo 
dominated nature of the present models. In other words, the 
dominant density distribution sets the natural scales. The 
easiest way to compute the DF is to change the integration 
variable from the total potential to the radius in the first of 
the identities in eq. ([2]). After normalization we obtain 



dF{v) 
dv 



v = v{Q), 



(27) 



where 



F(u) 



dg„ 



ds 



F,{») 



-*(«) 



and from eq. (f25 



1 



(28) 



< v < 00. (29) 



In the equations above, a negative sign appears in front of 
the integral (I28|) due to the monotonic decrease of the rel- 
ative potential with increasing radius; §„ is the normalized 
augmented density associated with eq. (|26jl . while s a = r a /ri 
is the normalized anisotropy radius; finally, the subscripts of 
the functions F refer to the isotropic and anisotropic parts 
of the DF, respectively. It is now evident how in halo domi- 
nated (n-7,71) models the external slope parameter appears 
in the density g n only, so that the DF of the n-7 component 
can be obtained by repeated differentiation with respect to 
£ of the simpler DF of the halo dominated (7,71) model. 
In the following, we will also discuss black hole dominated 
(n-7,BH) models: then, we will adopt as natural normaliza- 
tion scales the black hole mass Mbh and the scale-length r c 
of the n-7 component; moreover, < Q < 00, since in this 
case ^ = 1/s, so that in eq. (|29[1 Q = 

It is not difficult to show that, for generic values of 7, 
the integral (I28|l can be expressed as a combination of hy- 
pergeometric 2-F1 functions when the halo is a 71 = 0, or 1 
model, or a black hole. However, the resulting expressions 
are not more illuminating than the integral itself, so that 
we do not report them here. Instead, we prefer to show sim- 
ple asymptotic expansions relative to a couple of interesting 
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situations. The first concerns the behaviour of the DF of 
(n-7,BH) models for Q — > oo, for which the leading term of 
the expansion (normalized to /n/\/87I" 2 ) is 



f(Q) ~ 7 



7 ~2 



S7-3/2 



(30) 



thus, the slope 7 = 3/2 marks the different behaviour of the 
models. We note also that / < for 7 < 1/2 and the model 
is inconsistent (Tremaine et al. 1994); furthermore, we note 
how in the leading term above the function F a is producing 
no contribution. The other case of interest is the behaviour 
of the DF of isotropic n-7 models for Q — > 0. We recall 
that this case is very general, as it applies to self-gravitating 
n-7 models but also to n-7 models in generic halos of finite 
total mass, because in all systems *&t ~ GMt/s for s — * 00. 
In this case we found 



HQ) 



Q 



n-3/2 



(31) 



Fi- 



an expression (obviously) formally identical to eq. (|30[) 
nally, we notice that the asymptotic expansion for r — > 00 
of the velocity dispersion of isotropic n-7 models (both sin- 
gle component and embedded in dark matter halos of finite 
total mass) is of ~ l/( n + l) s i normalized to GMx/r c , so 
that for r — > 00 p/crf ~ (n + l) _3 / 2 s 3 / 2_n ; in accordance 
with the comment in Footnote 1. 



t<0 



f>0 

1 

WSC<0 




4.5 



5 

n 



5.5 



Figure 4. Upper limits for consistency on the concentration pa- 
rameter £ = r c /ri of the fully isotropic halo-dominated (n-0,1) 
model, as a function of the external density slope n. The solid 
line is the limit derived from the WSC (Section 3.2), while the 
dotted line is the true limit obtained from the DF (Section 4.3). 



holds, with 



4.1 The halo dominated (n-2,1) Model 

Following the preliminary discussion, here we derive the ex- 
plicit expression for the DF of the n-2 model with an arbi- 
trary degree of OM orbital anisotropy, in the gravitational 
field of a dominant Hernquist halo. After partial fraction de- 
composition of eq. (|28p . the isotropic and anisotropic com- 
ponents of the DF are 



and 



27r(n-4)! 



n — 3 (7 n — 4 



I Cjt3 + (jr 3 (_r2 — Or 2 



d £n-4 y £2 



(32) 



2yr(n-4)! d^™- 4 
Mv / TT^(n-2)!C n - 3 G„_i 



4vr(n-4)! 
respectively, where for k > 2 we define 



(33) 



G fc (£,1/) 



ds 



s-v(t+sy 



G° k ( V ) = G k (0, v). 



For ^ = 1 the function Gk assume a very simple form: 

r h ,A - v / ^ r ( fc - x ) 
k( ' >- r(fc-i/2)(i + I ,) fc - 1 ' 

while for £ 7^ 1 the formula 

ldGfctt.i/) _ (-l) fc - 1 d fc - 1 G 2 (^,i/) 



G fc+1 (4,^) 



fc! 



(34) 



(35) 



(36) 



+ 



1 + ;/ 



1 



arctan 



arctanh 



ll 


-€ 




f 1/' 


It 


- 1 



(37) 



(C > I)- 



Note that the second expression in eq. (|33[1 has been ob- 
tained from the relation (|36|) . From the derived formulae 
we determined the true anisotropy limit for the (n-2,1) mod- 
els. We found that /; > for all the explored values of n, 
i.e. the isotropic n-2 model in a dominant Hernquist halo 
is consistent independently of the halo scale-length, in ac- 
cordance with the result in point 1) of Section 3.2, obtained 
with the WSC. Therefore the set A+ in phase-space coin- 
cides with the whole accessible phase-space, and only the 
anisotropy limit s~ exists. In Fig. [S] (top panel) we show 
the isotropic (solid) and strongly anisotropic (i.e., near to 
the consistency limit) DFs for n = (4, 5, 6) and £ = 5. Note 
how the isotropic DFs are very similar, and monotonically 
decreasing. Instead, in the anisotropic cases the depression 
leading to inconsistency is apparent, clearly showing how 
phase-space inconsistency is always set outside the center. 
This feature is qualitatively similar to the others explored 
OM systems (e.g., see Fig. 2 in Ciotti & Lanzoni 1997, Figs. 2 
and 3 in C99, Fig. 3 in Ciotti, Morganti & de Zeeuw 2008); in 
addition, the systematic shift of the DF depression towards 
high Q values for increasing n is also apparent: an argument 
supporting this phenomenon is given in Section 2.1.1. We 
notice that the cuspy dips shown by the DFs for values of s a 
near the consistency limit could be the source of kinetic in- 
stabilities, whose investigation is of course well beyond the 
framework of this paper. In any case, it is almost certain 
that any system would develop radial orbit instability for 
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Figure 5. From top to bottom: the dimensionless DF for the halo 
dominated (n-2,1), (n-1,1), and (n-0,1) models, as a function of 
Q. The solid lines represent the isotropic DFs, while the dotted 
lines are the anisotropic DFs, with s a approaching the critical 
value for consistency. In all cases, £ = 5. 



Figure 6. Dimensionless DF for the black hole dominated (n- 
2,1) and (n-1,1) models, as a function of the relative energy Q, 
for n = 4, 5, and 6. Normalization scales are the mass of the 
black hole and the scale-length of the n-7 component. Solid lines 
represent the isotropic DFs, while dotted lines the anisotropic 
DFs, with s a approaching the critical value for consistency. 



values of the anisotropy radius larger than the consistency 
limit (e.g., see the N-body experiments discussed in Nipoti, 
Londrillo k, Ciotti 2002) : therefore, the critical DFs in Fig. 5 
describe equilibrium unstable systems. 

The trend of the anisotropy limit for different val- 
ues of n and increasing £ is shown in Fig. [3] (bottom panel) 
with the dashed lines. First, for fixed n and increasing £, 
s~ increases: in practice, for a given halo, a broader stel- 
lar density distribution is less and less able to sustain ra- 
dial anisotropy. Second, at fixed £ a steeper external density 
slope (i.e. larger n) corresponds to higher amount of admis- 
sible radial anisotropy, thus showing that not only the inner 
density slope, but also the external density profile is impor- 
tant for phase-space consistency. As expected, the true lim- 
itation on s~ is less stringent than the corresponding limit 
derived from the WSC, and the dashed line is everywhere 
below the solid line. In particular, while from the WSC the 
critical concentrations £ c in order to have a purely radial 
model are 1.5, 2, and 2.5 (for n — 4, 5, and 6), from the DF 
we obtain £ c ~ 2.8, 3.7, and 4.5. 

It is also possible to determine the DF of (n-2,BH) mod- 
els in presence of a dominant central black hole. The re- 
sulting formulae are identical to eqs. (|32ll - (|33| l. where the 
quantity yl + v at numerator is replaced by s/v and the G 
functions are replaced by the corresponding G. functions of 



same index, defined for k > 2 as 
s ds 



J V 



s-^(£ + s) 



(38) 



and 

GU»)^G. k{ O,„) = f 0^ . (39) 
For £ > the formula 

„ f t v ldG. h {t,v) {-If' 1 d^G^u) 



k d£ 



holds, where 



+ 



A-! 



arctanh 



d$ k 



(40) 



(41) 



Note that after performing the required differentiations 
one must set £ = 1. In Fig. [6] (upper panel) the obtained 
DFs are shown in the isotropic (solid lines) and strongly 
anisotropic (dotted lines) cases. Again, inconsistency due to 
OM anisotropy affects the DF at intermediate energies. The 
true anisotropy limits derived from the DF are given by solid 
squares in Fig. [2] Finally, note also the characteristic log-log 
linear trend of the DF for very high values of Q. This trend 
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is easily explained by considering the asymptotic expansion 
given in eq. ([30)l . 



4.2 The halo dominated (n-1,1) Model 

We now derive the explicit DF for the family of halo dom- 
inated (n-1,1) models, so that with this class of models we 
can study the effect of the external logarithmic slope when 
the two density profiles have the same inner slope. Of course, 
the (4-1,1) models are the halo dominated of (1,1) models 
discussed in C96. 

Integration of eq. (|28p reveals that the functions in- 
volved in the expression of the DF are the same as for the 
(n-2,1) models. In particular, 



Fi(u) = 



27r(n-4)! 
- 4 'G° 2 -G 2 



2G 3 _ 3G4.\ 



and 



FM = 



(-1)" wrr^ d"- 4 (2G 3 - 35G4) 



27r(n-4)! 



d^ n ~ 



(42) 



(43) 



In Fig. [5] (middle panel), the DF in the isotropic and strongly 
anisotropic cases is shown, for n = 4,5, 6. The behaviour is 
qualitatively similar to that of (n-2,1) models, i.e. while the 
isotropic systems display a monotonically decreasing DF, in 
the anisotropic case (when s~ approaches the consistency 
limit) the DF has a significant depression, that moves in- 
wards for increasing n. In Fig. [3] (middle row) we show the 
minimum value of the anisotropy radius normalized to r c as 
derived from the DF (dashed lines), and from the WSC. 
Again, in accordance with the preliminary study of Sec- 
tion 3.2, we found that f; > independently of the rela- 
tive concentration value £ and of the external slope n, so 
that only the limit s~ exists. All the trends exhibited by 
the (n-2,1) models are confirmed: in particular, sZ increases 
with £ for fixed n, while it decreases at fixed £ for increasing 
n. Note that the maximum radial anisotropy admissible for 
given £ and n is smaller than for (n-2,1) models, as expected 
from the shallower central profile of n-1 models. Finally, at 
variance with the (n-2,1) models (and in agreement with 
the expectations of the preliminary analysis) no purely ra- 
dial orbital configurations can be supported, independently 
of the value of the concentration £. 

As for (n-2,BH) models, also in the case of the black 
hole dominated (n-l,BH) models the functions F\ and F a 
are formally identical to those in eqs. (|42[) - (143p . with each 
G function replaced by the corresponding G* function, and 
with the substitution \/l + v — > ^Ji>; again, it is necessary 
to set £ = 1 after performing the required differentiations. 
While the true limits on the anisotropy radius are plotted in 
Fig.[2]with the solid squares, the obtained DFs are shown in 
the isotropic (solid) and strongly anisotropic (dotted) cases 
in the bottom panel of Fig. [6] not unexpectedly, the incon- 
sistency due to anisotropy manifests itself at intermediate 
energies. Finally, the behaviour of the DF for very high val- 
ues of Q is in accordance with the asymptotic expansion of 
the DF in eq. pO|l . 



4.3 The halo dominated (n-0,1) Model 

This case is expected to be the more complicated, because 
from C99 it is known that the (0,1) model presents a pecu- 
liar behaviour: in practice, there exists a range of £ where 
anisotropic models can be consistent when instead isotropic 
models are not, and this is due to the fact that /i becomes 
negative over the non empty set A-, so that two limits on 
s a must be considered. Note that the discussion of the WSC 
in the case of (n-0,1) models in Section 3.2 already showed 
a similar behaviour. 

As in the other cases, the integration of eq. (|28[) can be 
done by using the Gk functions only: 



(-l) n 3^/T+K 



n _ 3 d n-4 G; _ 



it (n - 



•4) 



n|3/n J +''C r 'G, 



7T 



(44) 



(-l) n 3liy/T+ZC 



2vr(n - 4)! 
" 4 (2£ 2 G 5 - 3£G 4 



+ G 3 



(45) 



where the second expression for Fi is obtained from the re- 
lation (|36[) of the Gfe functions. Motivated by the remarks 
above, we start to study the sign of fi, in order to deter- 
mine the condition for the existence of the sets A+ and A- . 
Indeed, we found that for £<5.233 (in the n = 4 case, see 
C99), £<6.192 (for n = 5), and £<7.166 (for n = 6), /; is 
positive everywhere over the accessible phase-space, a situ- 
ation which is graphically represented with the dotted lines 
in Fig. [4] It follows that for sufficiently concentrated models 
only the anisotropy limit s~ exists, and in Fig. [5] (bottom 
panel, solid lines) the isotropic DF is shown for (n-0,1) mod- 
els with £ = 5 and n = (4, 5, 6): note how, at variance with 
the other cases, the DF of (n-0,1) models decreases for in- 
creasing Q. The trend of s~ (£) is shown in Fig.[3](top panels, 
dashed lines): as for the models previously discussed, the Sa 
curve is similar (but displaced) with respect to that derived 
from the WSC. It is apparent how the s~ line lies above 
those of n-1 and n-2 models, as expected from the shallower 
density profile of n-0 models; moreover, we note how the 
minimum admissible values for decrease for increasing 
n, as in the other cases. In Fig. [5] the DFs for s a near the 
consistency limit s~ are shown with dotted lines for mod- 
els having £ = 5, and as in the previous 7 = 1 and 7 = 2 
models, inconsistency appears at intermediate values of Q. 

For values of £ larger than the critical values of £ c re- 
ported above, the set A- is not empty, because the isotropic 
component of the DF becomes negative at high relative ener- 
gies. In fact, at variance with the (n-2,1) and (n-1,1) models, 
the DF of the n-0 model converges to a finite limit (that can 
be easily calculated analytically from eq. [Jl]) for Q = 1; 
when increasing £ = ro/n, the value /(l) monotonically 
decreases, till it becomes negative for £ greater than £ c . A 
similar behavior was found in the numerical investigation 
of consistency of King (1972) and quasi-isothermal halos 
added to a de Vaucouleurs (1948) density distribution (see 
CP92). Thus, while the effects of OM radial anisotropy ap- 
pear to affect the DF at intermediate energies, the inconsis- 
tency of isotropic models embedded in a peaked halo seems 
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to characterize the DF at very high (relative) energies. Ac- 
cording to the discussion of Section 2, since the function / a 
is positive on A-, the upper bound s a can be determined: 
the appearence of this new value is manifested by the dotted 
lines in Fig. [3] with the vertical asymptote corresponding to 
£ c . Then, for £ greater than the critical limit, a range of con- 
centrations exists where an anisotropy radius r~ < r a < 
corresponds to consistent models; finally, for £ larger than 
the value where the dotted (sj) and dashed (s~) lines cross 
each other, s~ > s+, and the models become inconsistent. 
Summarizing, consistent models correspond to points placed 
above the s~ line and on the left of the s+ line. 

In the case of a central dominant black hole, the same 
comments as those for (n-2,BH) and (n-l,BH) models apply, 
and the functions F\ and -F a are obtained accordingly. Of 
course, the isotropic model is inconsistent, independently of 
the values of n, because the model with black hole would 
correspond to the formal limit £ — > oo in (n-0,1) models, 
which is already excluded by the previous discussion. What 
happens when considering the effect of anisotropy? Is there 
any compensating effect? A numerical inspection shows that 
the anisotropic part of the DF is negative on A-, so that 
not even orbital anisotropy can make the models consistent. 
Again, this result was expected from the study of (n-0,1) 
models, due to the crossing of the s~ and s a critical lines 
at finite values of f . 



5 DISCUSSION AND CONCLUSIONS 

In this paper, in a natural extension of previous investiga- 
tions (CP92, C96, C99), we focused on the importance of 
the external density logarithmic slope on the determination 
of the phase-space consistency of one and two-component 
stellar systems with OM orbital anisotropy. In particular, 
the considered (n-7, 71) models generalize (71,72) models, as 
the n-7 density component is characterized by logarithmic 
slope n instead of 4 outside the core radius, while in the 
inner regions the logarithmic density slope is 7. The main 
results can be summarized as follows: 

(i) It is shown that, for n integer, several structural and 
dynamical properties of the n-7 models can be obtained by 
differentiation, with respect to the scale-length, of the cor- 
responding formulae valid for the standard 7 models. For ex- 
ample, it is shown how to construct explicit solutions of the 
Jeans equations for the class of two-component anisotropic 
(711-71,712-72) models once the solution is known for (71,72) 
models. 

(ii) In one-component n-7 models, a lower limit for the 
anisotropy radius, so that smaller values certainly produce 
inconsistent models, is analytically derived following the 
technique introduced in CP92. It is found that for < 7 < 2 
this minimum anisotropy radius is strictly positive, but it 
decreases by increasing n, i.e. a larger amount of radial 
anisotropy can be supported by externally steeper density 
profiles. For 2 < 7 < 3 instead the necessary condition for 
consistency is satisfied Vs a > 0, independently of the value of 
7 and n. The minimum anisotropy radius so that larger val- 
ues correspond to consistent n-7 models is then determined 
by using the strong and weak sufficient conditions for model 
consistency derived in CP92. As for the necessary condition, 



we find that the minimum value of the anisotropy radius de- 
creases for increasing 7 and fixed n, and for increasing n and 
fixed 7. 

(iii) A similar analysis is then performed for two- 
component (n-7, 71) systems, in order to extend the study 
of C99. For simplicity we restrict to halo dominated models, 
i.e. we take into account only the gravitational field of the 
71 halo. In particular, we show that in the isotropic case, 
when 1 < 7 < 3 and < 71 < 7 (i.e., the halo density is 
centrally less peaked than the stellar component), the DF 
of the n-7 component is nowhere negative, independently of 
the mass and concentration of the 71 halo. This is true even 
when in the external parts the n-7 component is less peaked 
(3 < n < 4) than the 71 component. As a special application 
of this result, it follows that a black hole of any mass can 
be consistently added at the center of a globally isotropic 
n-7 model, when 1 < 7 < 3. Moreover, in the case of (n-0,1) 
models, where instead the halo is centrally steeper than the 
stellar component, the sufficient condition applied to the 
n-0 density distribution reveals the existence of an upper 
limit £ < (n + l)/2 of the relative concentration £ = ro/ri 
for model consistency in the isotropic case. Finally, in the 
case of anisotropic n-7 models with a dominant black hole 
at their center, we determined analytically a lower limit of 
the critical anisotropy radius for consistency as a function 
of n and 7; this value decreases for increasing 7 and/or n. 

(iv) The analytic expression for the DF of OM anisotropic 
halo dominated (n-7,1) models, and (n-7,BH) models with 
a dominant central BH is recovered in terms of elementary 
functions in the special cases when 7 = (0, 1, 2). It is found 
that, while for 7 = 1,2 the isotropic DF is positive inde- 
pendently of the model concentration, the 7 = component 
becomes inconsistent when the halo is sufficiently concen- 
trated, even in the isotropic case. In addition, the trend of 
the minimum value of the anisotropy radius as a function of 
the halo concentration confirms that steeper density profiles 
in the external region are consistent with a larger amount 
of radial anisotropy. 

In conclusion, the explored family of models provide a 
direct indication that the density slope in the external re- 
gions of stellar systems can be important in determining the 
admissible radial anisotropy, and this in addition to the well 
known relevance of the central density slope. A quantifica- 
tion of this argument is embodied in a necessary inequality 
for model consistency, that we derived by expressing a previ- 
ous finding (CP92) in terms of the logarithmic density slope. 
The obtained result holds for all consistent OM anisotropic 
systems, and relates the logarithmic density slope and the 
OM anisotropy indicator at each radius, in a way formally 
identical to the "cusp slope-central anisotropy" theorem by 
AE06, which is known to apply at the center of stellar sys- 
tems with generic anisotropy distribution, and everywhere 
in constant anisotropy systems. 
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APPENDIX A: CONSISTENCY 
REQUIREMENTS 

Al The NC for one— component n-7 models 

Simple algebra shows that the NC applied to anisotropic 
n-7 models can be written as 



2 - 7 - (n - 2)i 



s , < s < 00. 



(Al) 



7 + ns 

For 7 > 2 and n > 2 the r.h.s. is everywhere negative, so 
that inequality (Al) is trivially satisfied for s a > 0. Instead, 
for < 7 < 2 the parameter s a must be larger than the 
maximum of the expression of the r.h.s., which is reached at 



sm(h, 7) = 



n - (2n - 3)7 + y/(n - -y)[n + (4n - 9)7] 
2n(n - 2) 1 



(A2) 

that after substitution in eq. (Al) gives eq. (|21J1 . Of course, 
for n = 4 eq. (A2) in C99 is reobtained. 

A2 The WSC for halo— dominated (n-7 ,71) models 

Here we apply the WSC to the n-7 component of globally 
isotropic, halo-dominated (n-7,71) models, where 1 < 7 < 
3, < 71 < 7, and n > 3. Under the assumption of a dom- 
inant halo, eq. ((6]) reduces to investigate the positivity of 
an expression which factorizes in a strictly positive function 
and in an algebraic factor that after the natural substitu- 
tions 7 = 1 + e and 71 =7 — 61 (with < e < 2 and 
< ei < 7) becomes 

n(n - l)s 3 + n[n + 1 - e + ei + 2£e]s 2 (A3) 
+£{[(n + 1)(2 + e + ei) + eei - e] + fe 7 }s + f 2 7 (2 + 

whose positivity V(s,£) > is easily proved. Therefore the 
WSC is verified for any value of the concentration parameter 
£. Note that the variables £ and s adopted here correspond 
to l/£ and s/£ respectively in eq. (A13) of C99, due to the 
different normalization length. 

Being the WSC of the isotropic models satisfied, the 
existence of a critical anisotropy radius s~ for consistency 
depends on the sign of the anisotropic part of the WSC, as 
discussed in Section 2.2. In particular, the positivity of the 
anisotropic WSC for (n-2,1) models reads 



s(n-3)+n- l-2£> 0, < s < 00, 



(A4) 



and thus shows the existence of a critical value £ c = (n— 1) /2 
of the concentration parameter, marking an upper limit on 
the values of £ for which = 0, i.e. purely radial orbital 
distributions are allowed. 

The discussion of the WSC for globally isotropic n-0 
models in a dominant 7 = 1 halo is simple, reducing to the 
request that 



s(n - 1) + n + 1 - 2£ > 0, < s < 00, 



(A5) 



which is satisfied for £ < £ c = (n + l)/2. For n = 4 we 
reobtain the result of C99. 

In the case of central dominant BH, we assume in eq. (0 
Mt = Mbh, and from the previous discussion we restrict to 
n-7 models with 1 < 7 < 3. The WSC then requires that 

s 2 > g 2 (3 - 7) (7 - 2) + 2(n - 2)(3 - j)s - (n - 2)(n - 3)s 2 
n(n — l)s 2 + 2n(7 — l)s + 7(7 — 1) 



(A6) 

After the differentiation of the r.h.s. and the successive study 
of a quartic equation, it can be proved that for s > eq. (A6) 
admits only one maximum, located at sm = sivi(n,7) > 0. 
The general expression for sm(«, 7) is not reported here, but 
it was used to produce the long-dashed lines in Fig. [2] In 
any case, SM(n, 3) = 0, 



( 2 



s M (n,7) = < 



n — 3 ' 



/s 



2(n - 4) 



3n(n-l)(n-3) 3(n-l)(n-3) 
+ 4(n 3 +n 2 - 20n + 27) 



I 3(n-l)(n-3)^ 



where 

so = n [73n 4 - 660n 3 + 2262n 2 - 3484n + 2025 



(5n - 9)(13n 2 - 59n + 64)(n - l) 3 (n - 3) 3 



7 = 1, 

7 = 2, 
(A7) 

(A8) 



+ 9 



and for n = 4 and 7 = 1, 2, 3 the values given in Table 1 of 
C99 are reobtained. We finally report the reduced asymp- 
totic fourth-order equation needed to determine SM(n,7) in 
the limit n — + 00: remarkably, after the scaling s = y/n and 
the successive limit n — > 00, one obtains 

y 4 +2(2 7 -3)y 3 + 6(7-1) (7 ~2)j/ 2 (A9) 
+ 2(7 - 3)( 7 - 1)(2 7 - l)y + (7 - 3)( 7 - 2)( 7 - 1)7 = 0. 



APPENDIX B: THE VELOCITY DISPERSIONS 
AND THE VIRIAL QUANTITIES FOR 

(m-7i,n 2 -7 2 ) MODELS 

Here we present an easy way to express analytically the main 
dynamical quantities of the two-component OM anisotropic 
(ni-7i,n2-72) models for generic 71 and 72, and m and n2 in- 
teger > 4. The method is based on the evaluation of the sim- 
ilar quantities for the two-components (71,72) models, and 
then on repeated differentiation with respect to the scale- 
lengths. 

We start by considering the radial component a 2 of the 
velocity dispersion a 2 = <r 2 + a 2 of a system with a den- 
sity component p, which in the OM parameterization can 
be written as 



p(r)a T (r) 



where 



A(r)+r 2 J(r) 
r 2 + r\ 



(Bl) 



p(r)M T (r) 



p(r)M T (r)dr, I(r) =G n ''7" jr, 

(B2) 

(Binney & Mamon 1982), and Mx(r) is the total mass within 
r. Once of is known, the tangential velocity dispersion is 
obtained from eq. (|8]). 

As the method is general, we focus without loss of gen- 
erality on the component "1" of a two-component model, 
for which the two functions I and A are given by the sum 
Ji = In + J12 and Ai = An + An, due to the linearity 
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of mass in eqs. (B2), and where the meaning of the sub- 
script indices is apparent. Other quantities of interest in 
applications are the global energies entering the scalar virial 
theorem, and for the component 1 the scalar virial theorem 
reads 2Ki = |Wii| + |Wi 2 |, where Wn = -AttG f rpiMidr 
is the contribution due to the self-interaction, and W12 = 
—4nG J rp\Midr is the interaction energy with the "halo" 
(e.g., Ciotti 2000). 

The basic idea is to compute the integrals for the generic 
pair {p\,M2) relative to standard two-component (71,72) 
models with different scale-lengths r c i and r C 2, to perform 
the required differentiations as prescribed by eqs. (|18ft and 
Q19p . and finally (in the case of self-interactions) to set the 
two scale-lengths and the two slopes to the same value. From 
eqs. JTSJl and it follows 



dr, (B3) 



A l2 (r) = G 



r(m - 3)r(n 2 - 3) 

d" 1 " 4 d n ?- 4 f°° Pll {r,r cl )M 12 {r,r c2 ) 



dr 711 ' 4, dr 

III cl LLI : 



n 2 -4 
c2 



r c ir C 2 



"1-3 „"2-3 



(-1) 



r(m - 3)r(n 2 - 3) 

d™ 1 ' 4 d" 2 " 4 f 00 p 71 (r,r c i)Af 72 (r,r c2 ; 



Wia(r) = -4ttG 
d ni 



dr" 



-4 

c2 

,«l-3„ri! 
cl ' c2 



r c ir C 2r" 

3 /" i\rti+n 2 



dr, (B4) 



5 (-l) r 



r(ni - 3)r(n 2 - 3) 
1 cT 2 - 4 f°° p 71 (r,r c i)A/ 72 (r,r c2 ) 



r c ir C 2 



rdr;(B5) 



for generic values of 71 and 72 the integrals above involve 
hypergeometric 2-F1 functions, while for integer values of 71 
and 72 the resulting formulae can be expressed in terms of el- 
ementary functions (e.g., see C96 and C99 for (1,0) and (1,1) 
models; Ciotti et al. 1996 for (2,1) and (2,2) models). Even 
simpler expressions can be obtained without difficulty for 
to the halo dominated (n-0,1), (n-1,1), and (n-2,1) models, 
by adopting as normalization constants the physical scales 
Mi and n of the 7 = 1 dominant component, and by fixing 
ri2 = 4, r c2 = n, ni = n, £ = r c x/n, s = r/n, fx = M/Mi 
in eqs. (B3)-(B6), where M and r c i are the mass and the 
scale-lenght of the n-y component. As the integration is 
trivial but the results not particularly illuminating, we do 
not show the explicit formulae. 
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